Traffic jams induced by rare switching events in 
two-lane transport 



Tobias Reichenbach, Erwin Prey, and Thomas Franosch 

Arnold Sommcrfeld Center for Theoretical Physics (ASC) and Center for 
NanoScience (CeNS), Department of Physics, Ludwig-Maximilians-Universitat 
Miinchen, Theresienstrasse 37, D-80333 Miinchen, Germany 

E-mail: tobias . reichenbachSphysik . Imu . de 

Abstract. We investigate a model for driven exclusion processes where internal 
states are assigned to the particles. The latter account for diverse situations, ranging 
from spin states in spintronics to parallel lanes in intracellular or vehicular traffic. 
Introducing a coupling between the internal states by allowing particles to switch from 
one to another induces an intriguing polarization phenomenon. In a mesoscopic scaling, 
a rich stationary regime for the density profiles is discovered, with localized domain 
walls in the density profile of one of the internal states being feasible. We derive the 
shape of the density profiles as well as resulting phase diagrams analytically by a mean- 
field approximation and a continuum limit. Continuous as well as discontinuous lines 
of phase transition emerge, their intersections induce multicritical behavior. 
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1. Introduction 

Non-equilibrium critical phenomena arise in a broad variety of systems, including non- 
equilibrium growth models |1], percolation-like processes [2], kinetic Ising models [3], 
diffusion limited chemical reactions [4], and driven diffusive systems [5]. The latter 
provide models for transport processes ranging from biological systems, like the motion 
of ribosomes along a m-RNA chain or processive motors walking along cytoskeletal 
filaments [3, E], to vehicular traffic P, ^0\. In this work, we focus on the steady- 
state properties of such one-dimensional transport models, for which the Totally 
Asymmetric Simple Exclusion Process (TASEP) has emerged as a paradigm (for reviews 
see e.g. pT| [T2| [T3]). There, particles move unidirectionally from left to right on a one- 
dimensional lattice, interacting through on-site exclusion. The entrance/exit rates at 
the open left/right boundary control the system's behavior; tuning them, one encounters 
different non-equilibrium phases for the particle densities [r4 J. 

Intense theoretical research has been devoted to the classification of such non- 
equilibrium phenomena. For example, within the context of reaction-diffusion systems, 
there is strong evidence that phase transitions from an active to an absorbing state can 
be characterized in terms of only a few universality classes, the most important being 
the one of directed percolation (DP) |15]. To search for novel critical behavior, fruitful 
results have been obtained by coupling two reaction-diffusion systems [lEl |l7j, each 
undergoing the active to absorbing phase transition. Due to the coupling, the system 
exhibits a multicritical point with unusual critical behavior. 

We want to stress that already in equilibrium physics seminal insights have been 
gained by coupling identical systems. For instance, spin-ladders incorporate several 
Heisenberg spin chains [18|. There, quantum effects lead to a sensitive dependence on 
the chain number: for even ones a finite energy gap between the ground state and the 
lowest excitation emerges whereas gapless excitations dominate the low-temperature 
behavior if the number of spin chains is odd. 

In this work, we generalize the Totally Asymmetric Exclusion Process (TASEP) 
in a way that particles possess two internal states; we have recently published a short 
account of this work in Ref. [19]. Allowing particles to occasionally switch from one 
internal state to the other induces a coupling between the latter; indeed, the model may 
alternatively be regarded as two coupled TASEPs. When independent, each of them 
separately undergoes boundary-induced phase transitions [H] . The coupling is expected 
to induce novel phenomena, which are the subject of the present work. 

Exclusion is introduced by allowing multiple occupancy of lattice sites only if 
particles are in different internal states. Viewing the latter as spin-1/2 states, i.e spin- 
up (t) and spin-down (|), this directly translates into Pauli's exclusion principle; see 
Fig. [H Indeed, the exclusion process presented in this work may serve as a model for 
semiclassical transport in mesoscopic quantum systems [20], like hopping transport in 
chains of quantum dots in the presence of an applied field |21j. Our model incorporates 
the quantum nature of the particles through Pauli's exclusion principle, though phase 
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Figure 1. (Color online) Illustration of an exclusion model with two internal states, 
adopting the language of spin transport. Particles in states t (i) enter with rates 
{a^), move unidirectionally to the right within the lattice, may flip at rate oj and leave 
the system at rates (3^ {P^), always respecting Pauli's exclusion principle. 
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Figure 2. Illustration of the two-lane interpretation. We label the upper lane as lane 
I and the lower one as lane II. They possess individual entering rates, resp. a^^ as 
well as exiting rates, resp. 



coherence is ignored. A surprising analogy to a simple spintronics scheme, the Datta-Das 
spin field-effect transistor ^20j, holds. There, electrons move unidirectionally through a 
ferromagnetic metal or a semiconductor. The polarization of the electrons is controllable 
by a source for spin injection, a drain for spin extraction as well as a gate in the form 
of a tunable magnetic field that controls the strength of spin precession. In our model, 
this is mimicked by considering the spin-fiip rate as a control parameter. 

The model is potentially relevant within biological contexts, as well. In intracellular 
traffic [7[ [22] , molecular motors walking on parallel filaments may detach from one lane 
and attach on another, resulting in an effective switching between the lanes. In our 
model, identifying the two internal states with different lanes, one recovers a transport 
model on two lanes with simple site exclusion. In the same way, the system presented 
in this work serves as a highly simplified cartoon model of multi-lane highway traffic 
taking lane switching into account [9l [TO]. 

Significant insight into multi-lane traffic has been achieved (see Ref. [9l [10] and 
references therein). In particular, novel phases have been discovered in the case of 
indirect coupling, i.e. the velocity of the particles depends on the configuration on the 
neighboring lane [23[ [2H [25] . Recently, models have been presented that allow particles 
to switch between lanes, and the transport properties have in part been rationalized 
in terms of an effective single lane TASEP [26l [271 [2H]- There, the case of strong 
coupling has been investigated: the timescale of lane switching events is the same as of 
forward hopping. In our model, we explicitly want to ensure a competition between the 
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boundary processes and the switching between the internal states. We therefore employ 
a mesoscopic scaling, i.e. we consider the case where the switching events are rare as 
compared to forward hopping. This is the situation encountered in intracellular traffic 
^ where motors nearly exclusively remain on one lane and switch only very rarely. In 
the context of spin transport, it corresponds to the case where forward hopping occurs 
much faster than spin precession (weak external magnetic field). 

The outline of the present paper is the following. In Sec. [2] we introduce the model 
in the context of spin transport as well as two-lane traffic. Symmetries and currents 
are discussed, which play a key role in the following analysis. Section [3] describes in 
detail the mean-field approximation and the differential equations for the densities 
obtained therefrom through a continuum limit. The mesoscopic scaling is motivated 
and introduced, the details of the analytic solution for the spatial density profiles 
being condensed in Appendix A We obtain the generic form of the density profiles 
in Sec. HJ and compare our analytic results to stochastic simulations. We find that 
they agree excellently, suggesting the exactness of our analytic approach in the limit 
of large systems. As our main result, we encounter the polarization phenomenon, 
where the density profiles in the stationary non-equilibrium state exhibit localized 
'shocks'. Namely, the density of one spin state changes abruptly from low to high 
density. The origin of this phenomenon is rationalized in terms of singularities in 
coupled differential equations. We partition the full parameter space into three distinct 
regions, and observe a delocalization transition. The methods to calculate the phase 
boundaries analytically are developed simultaneously. Section [5] presents details on the 
stochastic simulations which we have carried out to corroborate our analytic approach. 
The central result of this work is then addressed in Sec. [6], where two-dimensional 
analytic phase diagrams are investigated. Our analytic approach identifies the phases 
where the polarization phenomenon occurs, as well as the continuous and discontinuous 
transitions that separate the phases. The nature of the transitions is explained by the 
injection/extraction limited current which is conserved along the track. As a second 
remarkable feature of the model, we uncover multi-critical points, i.e. points where 
two lines of phase boundaries intersect or the nature of a phase transition changes 
from discontinuous to a continuous one. Although multi-critical point are well-known 
in equilibrium statistical mechanics, a fundamental description for such a behavior for 
systems driven far from equilibrium still constitutes a major challenge. A brief summary 
and outlook concludes this work. 



2. The model 



In this section, we describe our model in terms of spin transport as well as two-lane 
traffic. Though we will preferentially use the language of spins in the subsequent 
sections, the two-lane interpretation is of no lesser interest, and straightforwardly 
obtained. Furthermore, we introduce two symmetries which are manifest on the level of 
the dynamical rules. 
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2.1. Dynamical rules 

We consider hopping transport on a one- dimensional lattice, composed of L sites, with 
open boundaries, see Fig. [H Particles possess internal states, which we restrict to 
two different kinds; adopting a spin notation, they are referred to as spin- up (t) and 
spin-down (|). They enter at the left boundary at rates resp. a^, and move 
unidirectionally from left to the right through the lattice. The timescale is fixed by 
putting the rate for these hopping events to unity. Within the bulk, particles may 
also flip their spin state, from spin-up to spin-down and back, at rate u. Finally, having 
reached the right boundary, particles may exit the system at rates (3^ resp. (3^, depending 
on their spin state. We allow all of these processes only under the constraint of Pauli's 
exclusion principle, meaning that every lattice site may at most be occupied by one 
particle of a given state. Spin-up and spin-down thus may simultaneously occupy the 
same site, however two particles with identical spin polarization cannot share a lattice 
site. In summary, our dynamical rules are the following: 

(i) at site i = 1 (left boundary), particles with spin-up (spin-down) may enter at rate 



(ii) at site i = L (right boundary), particles with spin- up (spin-down) leave the lattice 
at rate (3^ {(3^) 

(iii) particles may hop at unit rate from site i — 1 to the neighboring site i for 
i G {2, . . . , L}, i.e. within bulk 

(iv) within bulk, particles can flip their spin state with rate u, i.e. spin-up turns into 
spin-down and vice versa 

always respecting Pauli's exclusion principle. Processes (i) to (iii) constitute the Totally 
Asymmetric Simple Exclusion Processes (TASEP) for the two different states separately, 
while rule (iv) induces a coupling between them. Indeed, when the spin-flip rate uu 
vanishes, we recover the trivial situation of two independent TASEPs, while we will 
show that a proper treatment of u through a mesoscopic scaling induces nontrivial 
effects. 

2.2. Two-lane interpretation 

Having introduced our model in the language of semi-classical spin transport, where 
Pauli's exclusion principle is respected while phase coherence completely ignored, we 
now want to show that it also describes transport with site exclusion on two parallel 
lanes. As schematically drawn in Fig. [2], we consider two parallel lanes, each consisting 
of L sites, labeled as upper lane (I) and lower lane (II). They are identified with the 
internal states of the particles considered before: a particle with spin-up (spin-down) 
now corresponds to a particle on lane I (lane II). The processes (i) and (ii) describe 
entering of particles at lane I (II) at rate = (a^^ = a^) and exiting of lane I (II) at 
rate f3^ = {f3^^ = P^). Due to (iii), particles hop unidirectionally to the right on each 
individual lane; at rate u, they may switch from lane I to II and back. Pauli's exclusion 
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principle translates into simple site exclusion: all the above processes are allowed under 
the constraint of admitting at most one particle per site. Again, we clearly observe that 
it is process (iv) that couples two TASEPs, namely the ones on each individual lane, to 
each other. 

2.3. Symmetries 

Already on the level of the dynamical rules (i)-(iv) presented above, two symmetries are 
manifest that will prove helpful in the analysis of the system's behavior. We refer to 
the absence of particles with certain state as holes with the opposite respective state 0. 
Considering their motion, we observe that the dynamics of the holes is governed by 
the identical rules (i) to (iv), with "left" and "right" interchanged, i.e. with a discrete 
transformation of sites i ^ L — i a.s well as rates a^'^ ^ (3^'^ . The system thus exhibits 
a particle-hole symmetry. Even more intuitively, the two states behave qualitatively 
identical. Indeed, the system remains invariant upon changing spin-up to spin-down 
states and vice versa with a simultaneous interchange of ^ and (3^ ^ f3^, 
constituting a spin symmetry (in terms of the two-lane interpretation, it translates 
into a lane symmetry). 

When analyzing the system's behavior in the five-dimensional phase space, constituted 
of the entrance and exit rates a^"'-, f3'^'^ and u, these symmetries allow to connect 
different regions in phase space, and along the way to simplify the discussion. 

3. Mean-field equations, currents, and the continuum limit 

In this section, we shall make use of the dynamical rules introduced above to set up a 
quantitative description for the densities and currents in the system. Within a mean- 
field approximation, their time evolution is expressed through one-point functions only, 
namely the average occupations of a lattice site. Such mean-field approximations have 
been successfully applied to a variety of driven diffusive systems, see e.g. Ref. [12]. We 
focus on the properties of the non-equilibrium steady state, which results from bound- 
ary processes (entering and exiting events) as well as bulk ones (hopping and spin-flip 
events). Both types of processes compete if their time-scales are comparable; we ensure 
this condition by introducing a mesoscopic scaling for the spin flip rate uj. Our focus is 
on the hmit of large system sizes L, which is expected to single out distinct phases. To 
solve the resulting equations for the densities and currents, a continuum limit is then 
justified, and it suffices to consider the leading order in the small parameter, viz. the 
ratio of the lattice constant to system size. Such a mesoscopic scaling has been already 
successfully used in jSH [35] in the context of TASEP coupled to Langmuir dynamics. 

I The convention to flip the spin simuhaneously is natural in the language of solid-state physics. In 
the context of two-lane traffic, it appears more natural to consider vacancies moving on the same lane 
in the reverse direction. 
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3.1. Mean field approximation and currents 

Let n\{t) resp. be the fluctuating occupation number of site i for spin-up resp. 

spin-down state, i.e. n\'\t) = 1 if this site is occupied at time t by a particle with 
the specified spin state and n\'\t) = otherwise. Performing ensemble averages, the 
expected occupation, denoted by pj(t) and plit), is obtained. Within a mean-field 
approximation, higher order correlations between the occupation numbers are neglected, 
i.e. we impose the factorization approximation 

K{tXit))=plit)p^^{t); r,se{tl}. (1) 

Equations of motion for the densities can by obtained via balance equations: The 
time-change of the density at a certain site is related to appropriate currents. The 
spatially varying spin current jj{t) quantifies the rate at which particles of spin state 
t at site i — 1 hop to the neighboring site i. Within the mean-field approximation, Eq. 
([I]), the current is expressed in terms of densities as 

Jlit)=plim-Pi{t)], ^G{2,...,L}, (2) 

and similarly for the current The sum yields the total particle current Ji{t) = 

j]{t) + il{t). Due to the spin- flip process (iv), there also exists a leakage current jj^{t) 
from spin-up state to spin-down state. Within mean-field 

jI'it) = ujpim-pM, (3) 

and similarly for the leakage current ji\t) from spin-down to spin-up state. Now, for 
i G {2, . . . , L — 1} we can use balance equations to obtain the time evolution of the 
densities, 

^pIW =jl(t)-jl,it)+ji'(t)-jj'(t). (4) 

This constitutes an exact relation. Together with the mean field approximation for the 
currents, Eqs. ([2l[3]), one obtains a set of closed equations for the local densities 

^pIW = pUm - pM - pIW[1 - pIM + ^Plit) - ^pl(t) . (5) 

At the boundaries of the track, the corresponding expressions involve also the entrance 
and exit events, which are again treated in the spirit of a mean-field approach 

Apt (t)= a^l- pI (t)] - pI (t) [1 - pI (t)] + ujpi (t) - upl (t) , (6) 

j/lit) = pLiW[1 - Plit)] - + ^pi(t) - ujplit) . (7) 

Due to the spin symmetry, i.e. interchanging | and |, an analogous set of equations 
hods for the time evolution of the density of particles with spin-down state. 

In the stationary state, the densities p^ (t) do not depend on time t, such that the 
time derivatives in Eqs. dS])-© vanish. Therefrom, we immediately derive the spatial 
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conservation of the particle current: Indeed, summing Eq. (jlj) with the corresponding 
equation for the density of spin-down states yields 

J^ = Ji+l, te{2,...,L-l}, (8) 

such that the particle current does not depend on the spatial position i. Note that this 
does not apply to the individual spin currents, they do have a spatial dependence arising 
from the leakage currents. 

In a qualitative discussion, let us now anticipate the effects that arise from the 
non-conserved individual spin currents as well as from the conserved particle current. 
The latter has its analogy in TASEP, where the particle current is spatially conserved 
as well. It leads to two distinct regions in the parameter space: one where the current 
is determined by the left boundary, and the other where it is controlled by the right 
one. Both regions are connected by the discrete particle-hole symmetry. Thus, in 
general, discontinuous phase transitions arise when crossing the border from one region 
to the other. In our model, we will find similar behavior: the particle current is either 
determined by the left or by the right boundary. Again, both regions are connected by 
the discrete particle-hole symmetry, such that we expect discontinuous phase transitions 
at the border between both. Except for a small, particular region in the parameter space, 
this behavior is captured quantitatively by the mean-field approach and the subsequent 
analysis, which is further corroborated by stochastic simulations. The phenomena linked 
to the particular region will be presented elsewhere [29] . 

On the other hand, the non-conserved spin currents may be compared to the current 
in TASEP coupled to Langmuir kinetics; see Refs. [Ml [35]. Due to attachment and 
detachment processes, the in-lane current is only weakly conserved, allowing for a 
novel phenomena, namely phase separation into a low-density and a high-density region 
separated by a localized domain wall. The transitions to this phase are continuous 
considering the domain wall position as the order parameter. In our model, an 
analogous but even more intriguing phase will appear as well, with continuous transitions 
being possible. 

3.2. Mesoscopic scaling and the continuum limit 

3.2.1. Mesoscopic scaling. Phases and corresponding phase transitions are expected to 
emerge in the limit of large system size, L —>■ oo, which therefore constitutes the focus 
of this work. We expect interesting phase behavior to arise from the coupling of spin-up 
and spin-down states via spin-flip events, in addition to the entrance and exit processes. 
Clearly, if spin-flips occur on a fast time-scale, comparable to the hopping events, the 
spin degree of freedom is relaxed, such that the system's behavior is effectively the 
one of a TASEP. Previous work on related two-lane models [271 126] focused on the 
physics in that situation. In this work, we want to highlight the dynamical regime 
where coupling through spin- flips is present, however not sufficiently strong to relax the 
system's internal degree of freedom. In other words, we consider physical situations 
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where spin-flips occur on the same time-scale as the entrance/exit processes. Defining 
the gross spin-fiip rate Q = ujL yields a measure of how often a particle fiips its spin state 
while traversing the system. To ensure competition between spin-fiips with boundary 
processes, a mesoscopic scaling of the rate uj is employed by keeping fl fixed, of the same 
order as the entrance/exit rates, when the number of lattice sites becomes large L —>■ oo. 

3.2.2. Continuum limit and first order approximation. The total length of the lattice 
will be fixed to unity and one may define consistently the lattice constant e = 1/L. 
In the limit of large systems e — 0, a continuum limit is anticipated. We introduce 
continuous functions p\x) resp. p^{x) through p\xi) = p] resp. p^{xi) = pi at the 
discrete points Xj = ie. Expanding these to first order in the lattice constant, 

pTU) (x.iO = pT(i) (3;^ ± ,) = pTU) (3;^) ± ,d,ja) (3;^) , (9) 

the difference equations ©-(jT]) turn into differential equations. Observing that u = eQ 
is already of order e, we find that the zeroth order of Eq. ([5]) vanishes, and the first 
order in e yields 

[2p\x) - i]d^pHx) + np^{x) - np^x) = . (10) 

Similarly, the same manipulations for p^ yield 

[2p^{x) - l]d^p^ix) + np^x) - VLp^{x) = . (11) 
The expansion of Eqs. ([6]) and ([7]) in powers of e, yields in zeroth order 

p\0)= a\ p\l) = l-(3\ 

p^(0) = a^, p^(l) = 1 (12) 

which impose boundary conditions. Since two boundary conditions are enough to specify 
a solution of the coupled first order differential equations, the system is apparently 
over- determined. Of course, the full analytic solution, i.e. where all orders in e are 
incorporated, will be only piecewise given by the first-order approximation, Eqs. ( ITOi) - 
(fT2|) . Between these branches, the solution will depend on higher orders of e, therefore, 
these intermediate regions scale with order e and higher. They vanish in the limit of 
large systems, e — 0, yielding domain walls or boundary layers. 

Let us explain the latter terms. At the position of a domain wall, situated in bulk, the 
density changes its value discontinuously, from one of a low-density region to one of a 
high-density. Boundary layers are pinned to the boundaries of the system. There as well, 
the density changes discontinuously: from a value that is given by the corresponding 
boundary condition to that of a low- or high-density region which is imposed by the 
opposite boundary . 

3.2.3. Symmetries and currents revisited. In the following, we refiect important 
properties of the system, symmetries and currents, on the level of the first-order 
approximation, Eqs. (fT0l) - (fT2l) . The explicit solution of the latter can be found in 
Appendix A 
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The particle-hole symmetry, already inferred from the dynamical rules, now takes 
the form 

pTU)(a;) ^ 1 -p1^(^)(l -x), 

^ /3Ta) . (13) 

Interchanging | and | in the densities as well as the in- and outgoing rates yields the 
spin symmetry, 

T i 

^ /3^. (14) 

The individual spin currents as well as the particle current have been anticipated 
to provide further understanding of the system's behavior. In the continuum limit the 
zeroth order of the spin currents is found to be j^^^\x) = p^^-^\x)[l — p^^-^\x)], such that 
Eqs. ( JTOj) . ( ITTi) may be written in the form 

d^f = n[p^ - p^] , d,j^ = n[p^ -p^]. (15) 

The terms on the right-hand side, arising from the spin-flip process (iv), are seen to 
violate the spatial conservation of the spin currents. However, due to the mesoscopic 
scaling of the spin flip rate uj, the leakage currents between the spin states are only 
weak, see Eq. ([3]), locally tending to zero when e — > 0, such that the spin currents 
vary continuously in space. This finding imposes a condition for the transition from one 
branch of first-order solution to another, as described above: such a transition is only 
allowed when the corresponding spin currents are continuous at the transition point, 
thus singling out distinct positions for a possible transition. 

Finally, summing the two equations in Eq. ( |T5i) yields the spatial conservation of the 
particle current: d^J = 0. 

4. Partition of the parameter space and the generic density behavior 

The parameter space of our model, spanned by the five rates a^'^, P'^'^, and Q, is of 
high dimensionality. However, in this section, we show that it can be decomposed 
into only three basic distinct regions: the maximal- current region (MC) as well as the 
injection-limited (IN) and the extraction-limited one (EX). While trivial phase behavior 
occurs in the MC region, our focus is on the IN and EX region (connected by particle- 
hole symmetry), where a striking polarization phenomenon occurs. The generic phase 
behavior in these regions is derived, exhibiting this effect. 

4-1- Effective rates 

The entrance and exit rates as well as the carrying capacity of the bulk impose 
restrictions on the particle current. For example, the capacity of the bulk limits the 
individual spin currents j^^^'' to a maximal values of 1/4. The latter occurs at a density 
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Figure 3. Illustration of the effective rates. The right boundary is "open" , such 
that only the capacity of the bulk and the entrance rates limit the spin currents. The 
injection rate > ^ effectively acts as ^. The analytic predictions correspond to the 
solid lines, the results from stochastic simulations for L = 10000 are indicated by the 
wiggly line. With increasing spatial position, the densities approach a common value 
Pe- The parameters used are = 0.7, = 0.15, Q = 0.5. 



of 1/2, as seen from the previous result j^^-^-* = p^^-^^[l~ p^^^^)]. To illustrate the influence 
of the injection and extraction rates, we first consider an "open" right boundary i.e. 
fj'^ = fjl = I, Particles then leave the system unhindered, such that only the entrance 
rates may limit the particle current. Provided one of these rates, say a\ exceeds the 
value 1/2, the current of the corresponding state (j) is limited by the capacity of the 
bulk to a value of 1/4 in the vicinity of the left boundary. A boundary layer thus forms 
in the density profile of spin-up state at the left boundary, connecting the value of the 
injection rate to the value 1/2. Up to this boundary layer, the density profile p^(x) 
is identical to the one where takes a value of 1/2, c.f. Fig. [31 Similar reasoning holds 
for the extraction rates (3^^^\ They as well behave effectively as 1/2 when exceeding 
this value. To treat these findings properly, we introduce the effective rates 

1- 



T(i) 

a^g = mm 
iT(i) 



16a 



P;^^" = min[/3Ta),-J . (165) 

The system's bulk behavior will only depend on them, and, in particular, remain 
unaffected when a rate is varied at values exceeding 1/2. 



4-2. Injection-limited, extraction-limited, and maximal current region 

Equipped with these results, in the case of an "open" right boundary, the spin currents 
in the vicinity of the left boundary are given by = alg{l — alg) resp. j^ = alf^^l — alg), 
resulting in a particle current Jin imposed by the injection rates: Jin = «eff(-'- ~ 
alg) + alg{l — «cfr)- The analogous relations, with the injection and extraction rates 
interchanged, hold for the case of an "open" left boundary, = = 1. The particle 
current is then controlled by the right boundary: Jex = Plf^i^ — /^gg) + P^si^^ — Pis)- 
general, depending on which imposes the stronger restriction, either the left or the right 
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boundary limits the particle current: J < min( Jjn, Jex)- Indeed, J = min( Jjn, Jex) 
holds except for an anomalous situation, where the current is lower than this value [j. 
Depending on which of both cases applies, two complementary regions in phase space are 
distinguished: Jin < Jex is termed injection-limited region (IN), while Jin > Jex defines 
the extraction-limited region (EX). Since they are connected by discrete particle- hole 
symmetry, we expect discontinuous phase transitions across the border between both, 
to be referred as IN-EX boundary. 

Right at the IN-EX boundary, the system exhibits coexistence of low- and high- 
density phases, separated by domain walls. Interestingly, this phase coexistence emerges 
on both lanes (states), which may be seen as follows. Recall that a domain wall 
concatenates a region of low and another of high density. However, while the densities 
exhibit a discontinuity, the spin currents must be continuous. In other words, the spin 
currents, and therefore the particle currents, imposed by the left and right boundary 
must match each other. This yields the condition Jin = Jex, which is nothing but 
the relation describing the IN-EX boundary. Actually, what we have shown with this 
argument is that domain walls on both lanes (states) are at most feasible at the IN-EX 
boundary. However, it turns out that there, they do indeed form, and are delocalized. 
We refer to our forthcoming publication [29] for a detailed discussion of this phenomenon. 
Away from the IN-EX boundary, it follows that at most on one lane (state) a domain 
wall may appear. 

When both entrance rates a\ as well as both exit rates 13\ (3^ exceed the value 
1/2, the particle current is limited by neither boundary, but only through the carrying 
capacity of the bulk, restricting it to twice the maximal value 1/4 of the individual 
spin currents: J = 1/2. The latter situation therefore constitutes the maximal current 
region (MC). 



4-3. The generic state of the densities 

As we have seen in the previous section, particularly simple density profiles emerge in 
the MC region. There, up to boundary layers, the density profiles remain constant at 
a value 1/2 for each spin state. Another special region in parameter space is the IN- 
EX boundary, characterized by the simultaneous presence of domain walls in both spin 
states, as we discuss elsewhere [29] . 

Away from these regions, the generic situation for the density profiles is illustrated 
in Fig. m Here, we have considered parameters belonging to the IN region; the behavior 
in the EX region follows from particle-hole symmetry. A domain wall emerges for one 
spin state and a boundary layer for the other one. For specificity, we consider a domain 
wall for the spin-up state, the other situation is obtained from spin symmetry. The 

Tf I ) 

density profiles Pi close to the left boundary are given by the solution of the first-order 
differential equations (fTOj) . (fTTl) . obeying the left boundary conditions pj{x = 0) = 
and pI{x = 0) = a^g. For the density profiles pl^'^ in the vicinity of the right boundary, 

§ This situation arises in a certain neighborhood of the muhicritical points B, discussed in Sec. [5] 
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Figure 4. (Color online) The densities (a) and currents (b) in the IN region: generic 
state, exhibiting the polarization phenomenon. Results from stochastic simulations 
are shown as blue (L = 2000) resp. red (L = 10000) lines. They piecewise obey 
the first order approximations (black), grey lines indicate continuations of the latter 
into regions where the densities are no longer given by them. The parameters are 
= 0.4, = 0.2, (3^ = 0.2, (3^ ^ 0.45, and n = 0.5. 



we use the fact that the current in bulk is determined by the injection rates, J = Jin 
(which defines the IN region). Therefore, the densities satisfy right boundary conditions 
which are given by pj(a; = 1) = 1 — and pl{x = 1) is found from the conservation 
of the particle current: 

J = "Ifr(l - "Iff) + «Jff(l - «eff) 

= - Pis) +pUx = 1)[1 - pUx = 1)] . (17) 

At some point in bulk, the left and right solutions have to be concatenated by a 
domain wall for spin-up. To determine the position a;^ of this domain wall, we use the 
continuity of the spin currents; see Fig. |4](b). [J] This continuity condition singles out a 
distinct spatial position for the domain wall: Denote by p]{xw) the value of the density 
to the left of x^, and pl{xw) the value to the right. From = {1 — p^) together with 
pI{xw) 7^ pI{xw), we arrive at the condition 

pj(xu,) = 1 - pj(x^) (18) 

for the domain wall position. From the conservation of the particle current J, it 
follows that the density p^ is continuous at the position Xw 

II Indeed, though they are not spatially conserved, the mesoscopic scaling of the spin- flip rate lu was 
seen to cause a only slowly varying spatial dependence; in the continuum limit, the spin currents are 
continuous. 

% For TASEP-like transport the particle-hole symmetry restricts the density jump to this mirror 
relation. More general current-density relation are feasible [30l[3T], but are not expected to change the 
picture qualitatively. 
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When considering the internal states as actual spins, the appearance of a domain 
wall in the density profile of one of the spin states results in a spontaneous polarization 
phenomenon. Indeed, while both the density of spin-up and spin-down remain at 
comparable low values in the vicinity of the left boundary, this situation changes upon 
crossing the point Xw There, the density of spin-up jumps to a high value, while the 
density of spin-down remains at a low value, resulting in a polarization in this region. 

Comparing the generic phase behavior to the one of TASEP, we observe that the 
IN region can be seen as the analogue to the low-density region there: within both, 
a low-density phase accompanied by a boundary layer at the right boundary arises. 
Following these lines, the EX region has its analogue in the high-density region, while 
the MC region is straightforwardly generalized from the one of TASEP. Furthermore, 
the delocalization transition across the IN-EX boundary is similar to the appearance of 
a delocalized domain wall at the coexistence line in TASEP. 

4.4- Phases and phase boundaries 

In the generic situation of Fig. HI the density of spin-down is in a homogeneous low- 
density (LD) state, while for spin-up, a low-density and a high-density region coexist. 
We refer to the latter as the LD-HDin phase, as the phase separation arises within the 
IN region, to be contrasted from a LD-HDex phase which may arise within the EX 
region. Clearly, the LD-HDin phase is only present if the position of the domain 
wall lies within bulk. Tuning the system's parameter, it may leave the system through 
the left or right boundary, resulting in a homogeneous phase. Indeed, Xy^ = 1 marks 
the transition between the LD-HDin phase and the pure LD state, while at x^, = 
the density changes from the LD-HDin to a homogeneous high-density (HD) state. 
Regarding the domain wall position order parameter, these transitions are 

continuous. Implicit analytic expressions for these phase boundaries, derived in the 
following, are obtained from the first-order approximation, Eqs. (11 01) and ffTTj) . 

Spin symmetry yields the analogous situation with a domain wall appearing in 
the density profile of spin-down, while particle-hole symmetry maps it to the EX 
region, where a pure HD phase arises for one of the spins. Discontinuous transitions 
accompanied by delocalized domain walls appear at the submanifold of the IN-EX 
boundary (see [29] for a detailed discussion). 

The phase boundaries may be computed from the condition Xy^ = and x^, = 1 in 
the situation of Fig. HI Consider first the case of x^ = 0. There, the density profiles are 
fully given by the first-order approximation pr satisfying the boundary conditions at 
the right. The condition f|T8l) translates to 

pl(x = 0) = l-plix = 0) = l-al^ (19) 

which yields an additional constraint on the system's parameters. This defines the 
hyper-surface in the IN region where x^^, = occurs, and thus the phase boundary 
between the LD-HDin and the pure HD phase. 
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Similarly, if = 1, the densities follow the left solution (x), determined by the left 
boundary conditions, within the whole system. From Eq. ( ITSi) we obtain 

pj{x=l) = l-plix = l) = pls. (20) 

Again, the latter is a constraint on the parameters and defines the hyper-surface in the 
IN region where X uj — X IS found, being the phase boundary between the LD-HDin and 
the homogeneous LD phase. 

The conditions ( fT9l) . (l20l) yield implicit equations for the phase boundaries. The 
phase diagram is thus determined up to solving algebraic equations, which may be 
achieved numerically. Further insight concerning the phase boundaries is possible and 
may be obtained analytically, which we discuss next. 

First, we note that in the case of equal injection rates, = a^, the density profiles 
in the vicinity of the left boundary are constant. If in addition = = jS'^ < 1/2, 
we observe from Eq. fl20l) that a domain wall at = 1 emerges. Therefore, this set of 
parameters always lies on the phase boundary Xw = 1, independent of the value of fl. 

Second, we investigate the phase boundary determined by x^ = 0. Comparing with 
Fig. m we observe that the first-order approximation pj for the density of spin-up may 
reach the value | at a point which is denoted by xi: pl{xi) = |. This point corresponds 
to a branching point of the first-order solution. Increasing Q, the value of xi increases 

2 

as well. The domain wall in the density of spin- up can only emerge at a value x^, > xi. 
At most, which domain wall with infinitesimal small height arises. 

2 

For the phase boundary specified by x^ = 0, this implies that it only exists as long as 
xi < 0. The = corresponds to a domain wall of infinitesimal height, 

which is only feasible if a^g = ^. Now, for given rates = ^,a^,l3\ the condition 
xi = yields a critical rate fl*{a^,P^), depending on the rates a^,P^. The situation 
Xuj = can only emerge for rates Q < Q*{a^, (3^). Varying the rates a\a^ and f3\ the 
critical rate changes as well. In Appendix A we show that its largest value 



occurs ai = = 0. They yield the rate flc = = = 0), which is calculated 

to be 

fie = 1 + ^V21n(3 - 2v/2) ^ 0.38. (21) 

The critical Q*{a^,f3^) are lying in the interval between and Qc'- ^^*(a^,/5^) G [0,i7c7], 
and all values in this interval in fact occur. The rate Qc defines a scale in the spin- 
flip rate Q: For Q < Qq^ the phase boundary determined by x^ = exists, while 
disappearing for Q > Qc- 

Third, we study the form of the phase boundaries for large Q, meaning Q ^ Qc- 
In this case, the phase boundary specified by = is no longer present. Furthermore, 
it turns out that in this situation, the densities close to the left boundary quickly 
approximate a common value pe- The latter is found from conservation of the particle 
current: 2pe(l — Pe) = J- We now consider the implications for the phase boundary 
determined by Xw = 1. With pj{x = 1) = pe, Eq. ( l20l) turns into pe = PI^, yielding 

2/3jg(l - Pis) = alsil - als) + ai^il - ai^) . (22) 
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This condition specifies the phase boundary Xy^ = 1, asymptotically for large Q. It 
constitutes a simple quadratic equation in the in- and outgoing rates, independent of 
j3^, and contains the set = a^g = (3\ 

5. Stochastic simulations 

To confirm our analytic findings from the previous section, we have performed stochastic 
simulations. The dynamical rules (i)-(iv) described in Subsec. 12.11 were implemented 
using random sequential updating. In our simulations, we have performed averages 
over typically 10^ time steps, with 10 x L steps of updating between successive ones. 
Finite size scaling singles out the analytic solution in the limit of large system sizes, as 
exemplified in Figs. [3] and HI 

For all simulations, we have checked that the analytic predictions are recovered 
upon approaching the mesoscopic limit. We attribute the apparent exactness of our 
analytic approach in part to the exact current density relation in the steady state of the 
TASEP [32j. The additional coupling of the two TASEPs in our model is only weak: 
the local exchange between the two states vanishes in the limit of large system sizes. 
Correlations between them are washed out, and mean-field is recovered. 

The observed exactness of the analytic density profiles within the mesoscopic limit 
implies that our analytic approach yields exact phase diagrams as well. The latter are 
the subject of the subsequent section. 

6. Two-dimensional phase diagrams 

In this section, we discuss the phase behavior on two-dimensional cuts in the whole five- 
dimensional parameter space. Already the simplified situation of equal injection rates, 
= a^, yields interesting behavior. There as well as in the general case, we investigate 
the role of the spin-flip rate f2 by discussing the situation of small and large values of 

n. 

6.1. Equal injection rates 

For simplicity, we start our discussion of the phase diagram with equal injection rates, 
= a^. Then, the spin polarization phenomenon, depicted in Fig. HI becomes even 
more striking. Starting from equal densities at the left boundary, and hence zero 
polarization, spin polarization suddenly switches on at the domain wall position x^). 
The particular location of Xyj is not triggered by a cue on the track, but tuned through 
the model parameters. 

The phase transitions from LD to the LD-HDin arising in the IN region take a 
remarkably simple form. Their location is found from = 1, and is determined by 
Eq. ( l20l) (if phase coexistence arises for spin-up). Since p\x) = p^{x) = a = const, 
for a; < Xw, Eq. (l20l) turns into a = (3'^ . The latter transition line intersects the IN- EX 
boundary, given by Jin = Jex, at [3^ = [3^ = a, i.e. at the point where all entrance 
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Figure 5. Phase diagrams in the situation of equal entrance rates = = a and 
large f2. The phases of the densities of spin-up (spin-down) state are shown in (a) resp. 
(b) for a value = 0.3. At a multicritical point A, continuous lines (thin) intersect 
with a discontinuous line (bold), the IN-EX boundary. If > i, the maximal current 
phase appears for spin-up, see (c), as well for spin-down, drawn in (d). In the first 
situation, the switching rate is fi = 0.15, while f2 = 0.2 in the second. 



and exit rates coincide. At this multicritical point A, a continuous line intersects a 
discontinuous one. The same transition in the density of spin-down state is, from 
similar arguments, located at a = and also coincides with the IN-EX boundary 
in A. Neither the multicritical point A nor these phase boundaries depend on the 
magnitude of the gross spin flip rate VL. Therefore, qualitatively tuning the system's 
state is possible only upon changing the injection or extraction rates. The other phase 
transitions within the IN region, namely from the HD to the LD-HDin phase, are more 
involved. The analytic solution ( 1A.12I) . (1A.13I) has to be considered together with the 
condition f|T9l) for the transition. However, at the end of Subsec. 14.41 we have found 
that these transitions (determined by = 0) disappear for sufficiently large VL > flc- 

6.1.1. Large values ofQ. In the situation of large Q > flc, phase transitions arising 
from = in the IN region or from the analogue in the EX region do not emerge, as 
discussed at the end of Subsec. 14. 4[ We have drawn resulting phase diagrams in Fig. [5l 
showing the phase of spin-up (spin-down) in the left (right) panels, depending on a and 
Along the IN-EX boundary, being the same line (shown as bold) in the left and right 
panels, a delocalization transitions occur. At the multicritical point A, it is intersected 
by continuous lines emerging within the IN resp. the EX region. When f3^ > 1/2, a 
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Figure 6. Phase diagrams in the situation of Fig. [5] (a), (b), but with il decreased to 
a small value, = 0.05. Additional phase transitions emerge in the IN as well as the 
EX region, accompanied by multicritical points Bjj^ and B-^x- Caused by them, phase 
transitions do no longer appear across some parts of the IN-EX boundary, which is 
there shown as dashed line. 



maximal current (MC) phase emerges in the upper right quadrant, see Fig. O (c)-(d). 

To illustrate the system's phase behavior, let us consider what happens along a 
horizontal line in the phase diagrams (a) and (b), at a value > (3^ At such a line, 
for small values of a, both spin states are in low-density (LD) phases. Upon crossing 
a certain value of a, a domain wall enters aX = 1 in the spin-down density profile. 
Then, spin-down exhibits phase coexistence (LD-HDin), while spin-up remains in a LD 
phase. Further increasing a, the bold line is reached, where delocalized domain walls 
arise in both spin states. For larger values of a, a localized domain wall emerges for spin- 
up (implying a LD-HDex phase), and a pure HD phase for spin-down. If a is further 
increased, the domain wall in the spin-up density profile leaves the system through the 
left boundary (at Xy^ = 0), and pure HD phases remain for both spin states. 

While we have found the transitions within the IN region by simple expressions 
in the previous subsection, the ones emerging in the EX region are more complex and 
involve the full analytic solution flA.121) . flA.lSp . Their most notable feature is that the 



width of the corresponding coexistence phase decreases with increasing spin-flip rate 
Q, until it finally vanishes in the limit Q oo. This may be seen by considering the 
analogue of Eq. (122|) in the EX region, which describes the phase boundary as it is 
asymptotically approached when Q ^ oo: 

2aeff(l - «efr) = /3eff(l " PIs) + P^si^ ' /^eff) ! (23) 

it coincides with the IN-EX boundary. 



6.1.2. Small values ofVt. As discussed at the end of Subsec. 14. 4t when f2 < Vtc-, the 
appearance of additional phase transitions becomes possible. For example, within the IN 
region, the situation = may emerge. It describes the transition from the HD to the 
LD-HDiN phase; the analogue occurs in the EX region. In Fig. El we show resulting phase 
diagrams for the spin-up (left panel) and spin-down (right panel), resp.. The additional 
transition lines intersect the IN-EX boundary (bold) at additional multicritical points 
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i3iN and B^x- Also, they partly substitute the IN-EX boundary as a phase boundary: 
across some parts of the latter, phase transitions do not arise. This behavior reflects 
the decoupling of the two states for decreasing spin-flip rate Q. Indeed, for f2 — 0, the 
states become more and more decoupled, such that the IN-EX boundary, involving the 
combined entrance and exit rates of both states, loses its significance. 

6.1.3. Multicritical points. Although the shapes of most of the transition lines 
appearing in the phase diagrams shown in Fig. [6] are quite involved, they also exhibit 
simple behavior. Pairwise, namely one line from a transition in spin-up and another from 
a related transition in spin-down states, they intersect the IN-EX boundary in the same 
multicritical point. This intriguing phenomenon may be understood by considering the 
multicritical points: e.g., at A, the transition line from the LD to the LD-HDin phase in 
the density profile of spin-up intersects the IN-EX boundary, which implies that there 
we have a domain wall in the density profile of spin- up at the position = 1. However, 
being on the IN-EX boundary, the condition Jjn = Jex implies that in this situation a 
domain wall forms as well in the density of spin-down states, also located at a;^ = 1. 
Consequently, A also marks the point where the transition line specified by = 1 
for spin-down states intersects the IN-EX boundary. Due to the special situation of 
equal entrance rates, one more pair of lines intersects in this point. Similarly, at Bm, 
the transition line from the HD to LD-HDin phase in the density profile of spin-up 
intersects the IN-EX boundary, such that a domain wall forms in the density of spin-up 
at = 0. A gain, as Jjn = ^^ex holds on the IN-EX boundary, this implies the formation 
of a domain wall in the density of spin-down at x^; = 1, corresponding to the transition 
from the LD to the LD-HDex phase for spin-down in the EX region. 

6.2. The general case 

Having focused on the physically particularly enlightening case of equal entering rates in 
the previous subsection, we now turn to the general case. To illustrate our findings, we 
show phase diagrams depending on the injection and extraction rates for spin-up states, 
and Similar behavior as for equal entrance rates is observed. The multicritical 
point A now splits up into two distinct points and ^ex- 

6.2.1. Large values ofQ: Asymptotic results. Again, large Q prohibit the emergence 
of the phase transition from the HD to the LD-HDin phase in the IN region as well 
as from the LD to the LD-HDex phase within the EX region, see end of Subsec. 14.41 
In this paragraph, we consider phase diagrams which are approached asymptotically 
when fl ^ CO. Convergence is fast in Q, and the asymptotic phase boundaries yield an 
excellent approximation already for Q > 2Qc- 

The transition from LD to the LD-HDin phase in the IN region asymptotically 
takes the form of Eq. fl22|) . and the one from HD to the LD-HDex phase in the EX 
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Figure 7. Phase diagrams in the general situation: asymptotic results for large fl. 
Lines of continuous transitions (thin) within the IN resp. EX region intersect the 
delocalization transition line (bold) in multicritical points An<f resp. Aex- Both of 
these points appear in (a), (b) {a^ = 0.25, — 0.3) while only Aex is present in (c), 
(d) (a^ = 0.2,/?^ = 0.4) and Aw alone in (e), (f) {a^ = 0.05,/?^ = 0.15), yielding 
different topologies. 



region is obtained by particle-hole symmetry. All phase boundaries, including the IN- 
EX boundary, are thus given by simple quadratic expressions. 

Phase diagrams with different topologies that can emerge are exhibited in Fig. [7] 
and El As in the previous subsection, we show the phases of spin-up (spin-down) states 
on the left (right) panels. The phase boundaries between the LD and the LD-HDin 
phase in the IN region for spin-up as well as spin-down both intersect the IN-EX 
boundary in a multicritical point ^in, being located at (3^ = (3^ Similarly, the lines of 
continuous transitions within the EX region both coincide with the IN-EX boundary in 
a multicritical point ^ex, which is situated at = Note that the phase transitions 
emerging in the density profile of spin-down within the IN region do not depend on 
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Figure 8. Delocalization as well as maximal current phase (MC). When a-^ = (3^ < 
1/2 and a^,/?^ > 1/2 [upper right quadrant in (a), (b)], delocalized domain walls form 
in the density profiles of both spin states. If instead a^,/3-^ > 1/2, the maximal current 
phase emerges, see (c), (d). 

(3\ thus being horizontal hnes in the phase diagrams. Within the EX region they are 
independent of a\ yielding vertical lines. 

For a^,(3^ < 1/2, Fig. [7] shows different topologies of phase diagrams, which only 
depend on which of the multicritical points ^in, ^ex is present. If both appear, see 
Fig. (a), (b), the LD-HDin and the LD-HDex phase for spin- up are adjacent to each 
other, separated by the IN-EX boundary. Although in both phases localized domain 
walls emerge, their position changes discontinuously upon crossing the delocalization 
transition. E.g., starting within the LD-HDin phase, the domain wall delocalizes when 
approaching the IN-EX boundary, and, having crossed it, relocalizes again, but at a 
different position. 

When = (3^ < 1/2, a subtlety emerges, see Fig. [8] (a), (b). If both a'^ > 1/2 and 
> 1/2, i.e. in the upper right quadrant of the phase diagrams, these rates effectively 
act as 1/2, and the condition Jin = Jex for the IN-EX boundary is fulfilled in this whole 
region. Therefore, delocalized domain walls form on both lanes within this region, as is 
confirmed by our stochastic simulations [29] . 

The maximal current phase (MC) emerges when all rates exceed or equal the value 
1/2, corresponding to the upper left quadrant of the phase diagrams in Fig. [H] (c), (d). 




Figure 9. Phase diagrams in the general case and small values of fl. The nodal point 
A/jN remains unchanged when is varied. The appearance of the multicritical point 
Bij<j is accompanied by the non-occurrence of phase transitions across parts of the IN- 
EX boundary, then shown as dashed line. The multicritical point Ajj^x emerges in (a), 
(b), but not in the situation of (c)(, (d). Parameters are H. — 0.08, = 0.35, = 
0.45 in (a), (b) and Q = 0.2, = 0.15, (3^ = 0.4 in (c), (d). 



6.2.2. Small values ofVt. When Vt < Vie, the transitions from LD to LD-HDjn within 
the IN region as well as the analogue in the EX region are possible. As in the case 
of equal entering rates, the corresponding transition lines pairwise intersect the IN-EX 
boundary in multicritical points ;Bin and i^EX- As all transitions between phases of the 
spin-down density within the IN region are independent of (3\ the corresponding lines 
are simply horizontal; and within the EX region, their independence of a' implies that 
they yield vertical lines. The phase diagram for the density of spin-down is thus easily 
found from the IN-EX boundary given by Jin = </ex together with the locations of the 
multicritical points ^in, ^ex, and i^Ex- The latter follow from the intersection of 
phase transition lines for the density of spin-up, involving the whole analytic solution 
(IAI2D, (lAiall . with the IN-EX boundary. 

In Fig. [9] two interesting topologies that may arise are exemplified. Induced by the 
presence of the multicritical point Bm, phase transitions do not occur across all the 
IN-EX boundary, which is then only shown as dashed line. In Fig. [9] (a), (b), the points 
^Ex and i3iN are present. The LD-HDjn phase for spin- up intervenes the LD and the 
HD phase; the LD-HDex phase for spin- up is also present, though very tiny. In the 
phase diagram of spin-down, the LD-HDin phase intervenes the LD and the HD phase 
accompanied by continuous as well as discontinuous transitions. Again, the presence of 
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the multicritical points induces the topology; e.g. in Fig. M (c), (d), only Sin appears. 
For the discussion of the possible topologies, we encounter the restriction that Am and 
i3iN cannot occur together, as well as ^ex and Bex exclude each other (otherwise, the 
lines determined by a:^ = and = 1 would cross). 

We now discuss the influence of the spin-flip rate fl on the continuous transition 
lines for spin-up. In Subsec. I4.4l the manifold defined by = = a^g was found to be a 
sub-manifold of the phase boundary specified by = 1 in the IN region. Independent 
of Q, the point = = a^g, denoted by A/in, thus lies on the boundary between 
the LD and the LD-HDin phase (determined by = 1). For large Q, this boundary 
approaches the one given by Eq. (l22l) . 

Regarding the transition from the HD to the LD-HDin within the IN region (determined 
by Xu, = 0), Subsec. 14.41 revealed that for increasing Q it leaves the IN region at a critical 
transfer rate Q*{a^,P^). In the limit f2 — > 0, the densities p\x) and p^{x) approach 
constant values, and both the curve = 1 as = for spin-up in the IN region 
approach the line = for < |; The phase in the upper right quadrant in the 
phase diagram converges to the MC phase, such that in this limit, the case of two 
uncoupled TASEPs is recovered. 

7. Conclusions 

We have presented a detailed study of an exclusion process with internal states recently 
introduced in |19]. The Totally Asymmetric Exclusion Process (TASEP) has been 
generalized by assigning two internal states to the particles. Pauli's exclusion principle 
allows double occupation only for particles in different internal states. Occasional 
switches from one internal state to the other induce a coupling between the transport 
processes of the separate states. Such a dynamics encompasses diverse situations, 
ranging from vehicular traffic on multiple lanes to molecular motors walking on 
intracellular tracks and future spintronics devices. 

We have elaborated on the properties of the emerging non-equilibrium steady state 
focusing on density and current profiles. In a mesoscopic scaling of the switching rate 
between the internal states, nontrivial phenomena emerge. A localized domain wall in 
the density profile of one of the internal states induces a spontaneous polarization effect 
when viewing the internal states as spins. We provide an explanation based on the 
weakly conserved currents of the individual states and the current- density relations. A 
quantitative analytic description within a mean-field approximation and a continuum 
limit has been developed and solutions for the density and current profiles have been 
presented. A comparison with stochastic simulations revealed that our analytic approach 
becomes exact in the limit of large system sizes. We have attributed this remarkable 
finding to the exact current- density relation in the TASEP, supplemented by the locally 
weak coupling of the two TASEPs appearing in our model: u; — > in the limit of large 
system sizes. Local correlations between the two internal states are thus obliterated, as 
particles hop forward on a much faster timescale than they switch their internal state. 
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Furthermore, the parameter regions that allow for the formation of a localized 
domain wall have been considered. Analytic phase diagrams for various scenarios, in 
particular the case of equal entrance rates, have been derived. The phase diagrams 
have been found to exhibit a rich structure, with continuous as well as discontinuous 
non-equilibrium phase transitions. The discontinuous one originates in the conserved 
particle current, which is either limited by injection or extraction of particles. At the 
discontinuous transition between both regimes, delocalized domain walls emerge in the 
density profiles of both internal states. Multicritical points appear at the intersections 
of different transition lines organizing the topology of the phase diagrams. Two classes 
of multicritical points are identified, one of them arises only for sufficiently small gross 
spin-flip rate Q < Qc- The value Qc^ calculated analytically, provides a natural scale 
for the rate Q. 

It would be of interest to see which of the described phenomena qualitatively remain 
when generalizing the model to include more than two internal states. Indeed, within 
the context of molecular motors walking on microtubuli [7], between 12 and 14 parallel 
lanes are relevant. Also, the internal states might differ in the sense of different switching 
rates from one to another [28] and the built-in asymmetry may result in different 
phases. In the context of intracellular transport it appears worthwhile to investigate 
the consequences of a coupling to a bulk reservoir, c.f. [33l [Ml [35]; in particular, to 
study the interplay of domain wall formation induced by attachment and detachment 
processes as well as rare switching events. 
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Appendix A. The densities in the first order approximation and the critical 
value Qc 

In this Appendix, we give details on the derivation of the analytic solution of the 
mean-field approximation in the continuum limit to first order in e, i.e. the system 
of differential equations (fTOj) . (fTTI) . 
Summing them we find 



(A.l) 



such that 



[2p\xi)-l]' + [2p^{x,)-l]' 



J, 



(A.2) 
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constitutes a first integral. Remember that = p^{xi)[l — p^{xi)] + p^{xi)[l — p^{xi)], 
such that J is given by the total current: 

J = 2 - 4/°* . (A.3) 

This equation suggests the following parameterization: 

cos^ = J-5(2p^ - 1) 

sin^= J-5(2pi-i). (A.4) 
The derivative reads 



^/J je dp^ 

cos 0— = — — , 

2 ax ax 

which leads to the differential equation 

\fj sin Q cos Q— — 0(sin Q — cos 6) . 
dx 

This may be solved by a separation of variables: 



-X 



^(^^ sine cos ^ 



de. 



(A.5) 



(A.6) 



(A.7) 



^/J Jg(o) sin 61 -COS 61' 

To perform the integral, the substitution y — tan | is useful. We obtain the inverse 
function x — x{6): 



-G{y) 



+ 1. 



(A.8) 



j/=tan I 

Here we defined the function G{y) by 

and / is an constant of integration. 

To obtain the inverse functions x{p^) and x{p^), we have to express tan | by p^ 
resp. pK Recognize that tan| can be positive or negative. We therefore define 



-1 ifp^<i 
+1 if > I ' 

and analogously with | and I interchanged. Now 

e 



(A.IO) 



tan 



1 - J-i/2(2pt - 1) 
1 + J-V2(2pT _ 1) • 

The inverse functions x{p^) and x{p^) thus read: 



y=si 



y=s' 



l-J-l/2(2pT-l) 
l+j-l/2(2pT-l) 



l-J-l/2(2pi-l) 
l+j-l/2(2pi_i) 



(A.ll) 

(A.12) 
(A.13) 
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The constants of integration / and J = 2 — 4j*°* are determined by matching the bound- 
ary conditions. The inverse functions of Eqs. flA.121) . flA.13|) constitute the solution to 
Eqs. ( ITOl) and ( fTTi) . within the first-order approximation to the mean-field equations for 
the densities in the continuum limit. 



Next, we derive the result on Qq given at the end of Subsec. 14.41 Therefore, consider 
Fig. m We are interested in the point xi, and thus in the right branch of the spin-up 

2 

density profile. As the spin-down density is in the LD phase, i.e. it is smaller than i, we 
have = — 1 in the above solution for p^. Thus, y < in Eq. ( ]A.12I) . At the branching 
point of the analytic solution, i.e. the point xi, we have the density i, such that there 

2 ^ 

y = —1, implying G{y) = 0. Now, if this branching point lies on the right boundary, 
xi = 0, as it does for the critical fl*, this yields / = in Eq. flA.12p . On the other hand, 
the right branch satisfies the boundary condition on the right: p\x = 1) = 1 — f3\ 
Upon substitution into Eq. flA.121) . we obtain 

I = x{l - P^) = ^G{y) TV- (A.14) 



,/-l/2(l-2/3T) 
l+.J^l/2(l-2/3T) 



which for given a^, is an equation for Q*{a^, (3^). In particular, Q*{a^,P^) is 
monotonically increasing in G. Investigating G{y), it turns out that G{y) is in turn 
increasing in y. Since y is bounded from above by ?/ = 0, the maximal value for G{y) 
is provided by G{y = 0) = 1 + i\/2 In (3 - 2\/2). Next, we note that n*{a^,(3^) is 
an increasing function of y/j. With the constraint that = |, which is necessary 
for xi = 0, the largest y/j arises for = 0, i.e. y/j = 1. Combining both results, 
the maximal value for the critical rates fl*{a^,(3'^) occurs for = and y = 0. Both 
conditions together yield 

f]c = l + ^V21n(3-2v^). (A.15) 

Finally, we note that = and y = implies (3^ = 0, such that Qc arises if a-'- = 
and (3^ = 0. 
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